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Abstract 

We present a new method for proving non-holonomicity of sequences, which 
is based on results about the number of zeros of elementary and of analytic 
functions. Our approach is applicable to sequences that are defined as the values 
of an elementary function at positive integral arguments. We generalize several 
recent results; e.g., non-holonomicity of the logarithmic sequence is extended 
to rational functions involving log n. Moreover, we show that the sequence that 
arises from evaluating the Riemann zeta function at odd integers is not holonomic. 

Mathematics Subject Classifications: 11B37, 11B83 

1 Introduction 

A sequence of complex numbers (a n ) n >i is called holonomic (P-recursive) if there exist 
polynomials Pk(x) € C[x], k = 0, 1, . . . , d, which are not all zero, such that 

p (n)a n + p 1 (n)a n+1 H \-p d (n)a n+d = (1) 

holds for every integer n > 1 (equivalently, for every integer n > no). Their continuous 
counterpart are holonomic (D-finite) functions g(z), which satisfy linear differential 
equations 

Po{z)g(z)+ Pl {z)g'(z) + ■ ■ ■ + p d {z)g {d \z) = 

with polynomial coefficients. For more information on holonomic sequences and func- 
tions, especially in combinatorial enumeration, see Stanley |21| . For the sake of clarity, 
in this paper we will use the word holonomic for sequences and the word D-finite for 
functions. Sequences that satisfy Q with d = 1 are called hypergeometric. An im- 
portant property for our purposes is that the sum and the product of two holonomic 
sequences are holonomic. 

The role of the holonomic sequences in the set of all complex sequences is reminis- 
cent of the role of the algebraic numbers in the set of all complex numbers. Just like 

*Simon Fraser University, Vancouver, jpb at sfu.ca 

1 Vienna University of Technology, sgerhold at f am. tuwien. ac . at 

t Charles University, Prague, klazar at kam.mff.cuni.cz 

^Institute de Matematicas UNAM, Campus Morelia. fluca at matmor . unam . mx 



1 



complex numbers are usually not algebraic, unless they are by design, a sequence that 
is not obviously holonomic is usually not holonomic. In both situations it can be a 
challenge, though, to come up with a proof. Flajolet, Gerhold, and Salvy ^ give an 
exhaustive survey of known non-holonomicity results. In the present article we shall 
be interested in the following 

Question. If / : — > C is a "nice" function that is defined on a domain O C C 
containing all natural numbers N = {1, 2, 3, . . . }, when is it the case that the sequence 
of values of / on N, (/(n))„>i, is non-holonomic? 

We seek general results and criteria yielding non-holonomicity of (/(«)) n >i for / from 
large classes of functions J 7 , generalizing the results from Gerhold [Sj and some of the 
results of Flajolet, Gerhold, and Salvy 0]. In the first paper, (logn)„>i was proved 
conditionally to be non-holonomic, and in the second paper an unconditional proof 
was given by an asymptotic machinery. A simple proof was given also by Klazar 
Here, we generalize the method from the latter paper and arrive at non-holonomicity 
results for many sequences for which closed-form expressions are available. 

In the following section, we explain our proving method. Section [3] is devoted to 
straightforward applications of our approach. In Section 0] we develop our arguments 
further in several directions to establish the non-holonomicity of some sequences in- 
volving the exponential and factorial functions. We proceed by presenting several 
unrelated results, e.g. about algebraic sequences and interlacement sequences, in Sec- 
tion Finally, we show in Section that the sequence arising from evaluating the 
Riemann zeta function at odd integers is not holonomic. The method we use there is 
not related to our main approach. 

2 The Proving Method 

Before presenting our approach we briefly discuss the method of Flajolet, Gerhold, and 
Salvy They use the equivalence saying that (o„)„>i is holonomic iff the ordinary 
generating function g(z) = J2 n a nZ n is D-finite [23- The asymptotic behavior of a 
holonomic function near a singularity is constrained by a structure theorem 4, The- 
orem 2], and Abelian theorems transfer the asymptotic behavior of a n as n — > oo to 
the asymptotic behavior of g(z) near the singularity Theorem 3]. Thus, if (a n ) n >i 
(or any transform of it obtained by holonomicity-preserving tranformations) has an 
asymptotic behaviour that is transferred to "forbidden" asymptotics of the generat- 
ing function, then (a„)„>i is not holonomic. While requiring an initial asymptotic 
analysis of (a n )n>ii which is a drawback, the real strength of this method is that this 
asymptotic behaviour is fully sufficient information, and a closed-form representation 
a n = f(n) is not needed to prove non-holonomicity. This has been demonstrated, e.g., 
on the sequence of primes @|- In contrast, our method relies heavily on the explicit 
representation a n = f[n). 

We associate with a function f[z) and a [d + l)-tuple of complex polynomials 
Po(z), . . . ,pd(z) the function 

d 

F(z) = F(z; f, Po , ...,Pd):= + *)■ 

i=0 
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If (f(n)) n >i is holonomic, then, for some not-all-zero polynomials po, . . . ,pd, the func- 
tion F(z) vanishes at z = 1, 2, If / is real-valued in the positive reals (or at least in 

N), then (/(n)) ra >i is a sequence of real numbers, and by an algebraic argument [Hl ll^j 
the polynomials in Q can be assumed without loss of generality to have real coeffi- 
cients. This is important for the formulation of the following two properties for classes 
of complex functions T: If all / G T are real-valued for real arguments, we assume 
that the polynomials po(z), . . . ,Pd(z) have real coefficients. 

Property (A) If / G T is not identically zero and not all (real or complex, see above) 
polynomials po(z), . . . ,pd{z) are zero, then the function F(z; f,po, ■ ■ ■ ,Pd) is not 
identically zero. 

Property (B) If F(z) = F(z\ f,po, ■ ■ ■ ,Pd), f £ J 7 , vanishes at all z G N, then F(z) 
is identically zero. 

It follows that any non-identically-zero function / from a class T with properties (A) 
and (B) produces a non-holonomic sequence (/(n))„>i. (For some of the classes T we 
are going to discuss, 1 has to be replaced by some other positive integer no.) 

Let us have a look at some (classes of) functions having property (A) and then at 
functions having property (B). Ultimately we want to have large classes of functions 
having simultaneously properties (A) and (B). 

Condition (Al). Meromorphic functions and singularities. The class 
of functions / which are meromorphic for $t(z) > but have no extension that is 
meromorphic at z = has property (A). Indeed, if po is a nonzero polynomial and 
F(z; f,p , ...,p d ) = 0, then 

d 

f(z)^^p (z)- 1 Y / M^f^ + ^) (2) 

for 5R(z) > 0, which is impossible, because the right hand side is meromorphic at z — 
while the left hand side is not. This criterion can be generalized in an obvious way by 
replacing with any other point zq £ C. Examples of such functions are log z, log log z, 
z a with a € C\Z, exp(- v /z), exp(-\/log z), sin(l/z), z z , arctanz, etc. Condition (Al) 
will be our main tool for establishing property (A). Note how the assumption that our 
sequences a n are defined via functions f(z) enables us to work with the asymptotic 
behaviour of f(z) at points other than infinity. 

Condition (A2). Logarithmic derivatives and shifts. Here we assume that 
f'(z)/f(z) £ C(z). Notice that this implies that / is of the form 

m 

f(z) = exp(r(z)) J[(z - a*) 04 , 

i=l 

where r(z) is a rational function, and cti, Ci G C. 

Proposition 1. The class of functions f(z) that are analytic in a neighborhood of 
[l,oo[ and satisfy f'(z)/f(z) G C(z) but f(z + k)/f(z) G" C(z), k G N, has property 
(A)- 
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Proof. Suppose, on the contrary, that we have a nontrivial relation F[z\ f,po, ■ ■ ■ , pa) = 
0, with minimal d. Upon differentiating, this implies 

o = f'(z) = J2 (p'^+p^j^rf) ^ z + *)• 

We now eliminate the summand i — -- from the equations F{z) = and F'{z) = 0. 
The coefficients b%(z) in 

J2 bi(z)f(z + i) := F'(z)p (z) - F{z) (p' (z) + p (z)^j g C(z) 
are in C(z), hence they vanish by the minimality of d. For i = d this reads as 
0=Po(z) lp d (z)+p d (z) yr^-—K I - [Po{ z )+Po{z)j^j J Pd(z). 

Dividing by po(z)pd(z) and integrating, we obtain 

f(z + d)/f(z) = Cp (z)/p d (z)eC(z) 

for some constant C, a contradiction. Note that we have po(z)p ( i(z) by the 
minimality of d. □ 

Notice that the field C(z) in the proposition can be replaced by an arbitrary 
field K D C(z) of functions meromorphic at infinity, provided that it is closed un- 
der differentiation and shift. For instance, this generalization can be applied to 
f{z) = expexp(l/z) with K = C(z, e 1 ^ ,e 1 ^ z+1 \ . . .). 

We remark that we do not know of examples where condition (A2) works and 
condition (Al) fails. However, Proposition ^ can be viewed as a first step to address 
the question which D-finite functions f(z) define holonomic sequences (/(w)) n >i- The 
interplay of shift and derivative in the proof of Proposition ^ shows that this problem 
is not as contrived as it might seem. 

Condition (A3). Growth conditions. If f(z) is such that \f(z + l)/f(z)\ grows 
faster than polynomially as \z\ — > oo then J2J cannot hold identically either. We will 
not use this observation in what follows, since it is well-known that sequences that 
grow faster than a power of n! cannot be holonomic. 

We now present some sufficient conditions for property (B). 

Condition (Bl). Rational functions. The class of rational functions has prop- 
erty (B). This is trivial but by itself not very useful, since rational functions do not 
have property (A). Besides, rational functions clearly define holonomic sequences. 

Condition (B2). Differentiation and Rolle's theorem. If F(z) is a non- 
identically-zero real smooth function and has infinitely many real positive zeros, then 
by Rolle's theorem so have F'(z) and all higher derivatives of F(z). Thus, if for 
sufficiently large k we get a function F^ (z) with only finitely many real positive zeros 
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(e.g., a rational function, cf. condition (Bl)), then F(z) has only finitely many zeros 
as well. This argument shows that, for example, the class of functions that lie in 
ffi.[x, logx] and are defined in [l,oo[ has property (B). 

Condition (B3). Meromorphicity at infinity. The class of functions f(z) 
meromorphic at oo has property (B). Indeed, the function F{z) = F(z; f,po, ■ ■ ■ ,Pd) 
is then meromorphic at oo as well and cannot have infinitely many zeros in N, unless 
it is identically zero. This applies, e.g., to sin(l/z), exp(l/z), etc. 

Condition (B4). Classical results about zeros of analytic functions. Many 
authors have investigated the distribution of zeros of analytic functions. We just recall 
the following classical theorem |15|. Its gist is that the function sin7rz is the "smallest" 
function analytic in the right half-plane that vanishes at 0, 1, 2, ... . 

Theorem 2 (Carlson). Let g(z) be a function analytic in 5R(z) > that satisfies 
the growth conditions g(z) = 0(exp(a|z|)) and g(±ix) = 0(exp((7r — s)x)), x > 0, 
for some positive constants a and e. If g(z) vanishes at z = 0, 1, 2, ... , then g(z) is 
identically zero. 

Corollary 3. The class of functions f(z) such that g(z) = f(z + 1) satisfies the 
conditions of Theorem^ has property (B). 

Condition (B5). Zeros of real elementary functions results of Kho- 
vanskii. In his book Fewnomials ^UJi Khovanskii proves rather general results on 
the number of simultaneous zeros of multivariate real elementary functions. We will 
use his results only for univariate functions. Adopting Khovanskh's definition, we 
call a univariate function g{x) elementary ; if it can be expressed as a composition of 
multivariate rational functions and the functions exp(x), log(x), sin(x), cos(x), tana;, 
arcsin(a;), arccos(a;), and arctan(x). The domain of definition of g{x) must be such 
that denominators do not vanish, arguments of log are positive, arguments of arcsin 
and arcos lie in [—1, 1], arguments of tan in [— tt/2, 7r/2], and arguments of sin and cos 
are bounded. The latter assumption is crucial to recover the finiteness of the zero set. 
Then a special case of one of Khovanskh's main results ^3 §1.6] says: 

Theorem 4 (Khovanskii). An elementary function g{x) has only finitely many sim- 
ple zeros in its domain of definition. 

Though apparently only a technical nuisance, the restriction to simple zeros makes 
this result somewhat unwieldy for our purposes. To overcome this difficulty, we impose 
the additional assumption that f(x) be a quotient of D-finite functions. Then F(x) is 
also a quotient of D-finite functions, and it is easy to see that the order of the zeros of 
F(x) is therefore bounded. Thus, by repeated differentiation, we arrive at a function 
with only simple zeros, to which Theorem 01 can be applied. 

Corollary 5. The class of elementary functions that are defined on [l,oo[ and are a 
quotient of D- finite functions has property (B). 

Let us briefly compare this to Carlson's theorem. Apart from the absence of a 
growth condition, Khovanskh's result does not require analyticity in a right half-plane. 
Thus, for instance, we need not worry about zeros of the denominator of f{x) outside 
the positive reals. On the other hand, Khovanskh's result is about real functions only, 
and his restriction to simple zeros induced us to include the additional Z3-finiteness 
requirement. 
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3 First Results 



Putting together the results from the preceding section immediately yields several 
non-holonomicity proofs without any additional work. 

Theorem 6. Let f(z) be a function from M.(z, logz, arctanz) that is analytic in an 
open set containing the positive reals. Then (f(n)) n >i is a holonomic sequence if and 
only if f{z) is a rational function. 

Proof. The result is an immediate consequence of conditions (Al) and (B2) (Corol- 
lary0. Indeed, if f(z) depends on logz, then F(z) = F(z;f,p , . . . ,pd) has a loga- 
rithmic singularity at z = 0, and if it depends on arctanz, then F(z) has a logarithmic 
singularity at z = i. □ 

We have excluded the entire functions cxp and sin from the preceding theorem, 
since condition (Al) does not work for them in general. In Section0|we will completely 
describe the rational functions of n and e™ that define holonomic sequences. As long 
as there is some log or arctan present, the argument from the proof of TheorcmElgoes 
through: For instance, 

_ logn + sin^j 
a ' n ~ e™ + n 2 + 1 

does not define a holonomic sequence for this reason. As for compositions, functions 
like log log z are not allowed in Theorem|Hl since Infinite functions are not closed under 
composition, violating our Z?-finitcncss assumption that takes care of multiple zeros of 
elementary functions. (The sequence (loglogn)„>2 can be treated by the asymptotic 
method 0] for showing non-holonomicity.) Still, D-finite functions are closed under 
composition with algebraic functions. For example, 

ft \ &Z2 + 1 

f( z ) = — 7 

e 2 — 1 

is a quotient of D-finite functions, hence the assumptions of Corollary are satisfied. 
The complex poles of f(z) are at z = ±^/2kire 17T / 4 , k 6 Z, which shows that con- 
dition (Al) can be applied, establishing the non-holonomicity of the corresponding 
sequence. Observe that Carlson's theorem is not (directly) applicable here, since f(z) 
is not analytic in a right half-plane. Anyhow, let us now recapitulate what we can 
infer from Corollary 

Theorem 7. Let f{z) be a function analytic for 5ft(z) > that does not have an 
extension that is meromorphic at z = 0. If f(z+ 1) satisfies the growth conditions in 
Theorem^ then the sequence (f(n)) n >i is not holonomic. 

This gives another proof that logn is not holonomic, and also establishes non- 
holonomicity of e^™, e 1 /™, expexp(l/n), and n a with a S C \ Z. Note that for the 
latter three we could have used condition (A2) instead of (Al) to establish that F 
cannot vanish identically. Also, Khovanskh's theorem can be used instead of Carlson's 
in all five cases. Observe that, to our knowledge, this is the first published proof that 
e^/™ and expexp(l/n) are not holonomic. 

We have seen by now that our method works on several sequences without any 
effort. Take for example the sequence e 1 /™. Proving its non-holonomicity by asymp- 
totics requires a somewhat involved asymptotic analysis [J] of the generating function 
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Y^ n >i e 1 ^ 2 ™- But longer proofs can of course have merits, too; this analysis has stim- 
ulated comprehensive further asymptotic investigations jS] of analytic functions of the 
form X!n>i °™ w hich are of interest in their own right. 

4 Factorial and Exponential Sequences 

We now present some less obvious applications of the tools we have collected in Sec- 
tion [3 We start with some rather concrete sequences and give a few more general 
results, most notably about algebraic sequences, in the following section. Let us begin 
with the sequence nl a , aeC\l It grows too fast to apply Carlson's theorem. Also, 
the function T(z + l) a is not a quotient of D-finite functions, hence Corollary [3] is not 
applicable. We now show that this sequence can be knocked out by Carlson's theorem 
after a little bit of massage on the alleged recurrence. 

Theorem 8. Let u\,...,u s be distinct complex numbers, and let a±, ... , a s be complex 
numbers. Then the sequence T(n — u\) ai . . .T(n — u s ) Qs is holonomic if and only if 
otl, . . . ,a s are integers. 

In particular, powers of hypergeometric sequences with non-integral exponent are 
not holonomic, unless the exponent trivially cancels, such as in (n! 2 ) 1 / 2 . This gener- 
alizes a known result [SJ Theorem 1]. 

Proof. Suppose that the sequence is holonomic. Since hypergeometric sequences are 
holonomic, we may assume w.l.o.g. that no on is an integer. The sequence is defined for 
n > uq := [a\ + 1, where a :— max^ 3?(-Ui). Suppose that ^2t=oP k ( z )f( z + k) vanishes 
for z = 7101^0 + 1) •• •) where the Pk(z) are polynomials, po(z) is not identically zero, 
and 

f(z) :=r(z- Ul ) ai ...T(z~u s ) a °. 
By the recurrence of F(z) and the fact that T(z) has no zeros, this implies 

d s 4 

^Pk(z) J| fO - =0, z = n ,n + 1, . . . , (3) 

k=Q i=l 

where z k denotes the rising factorial 

z 1 = z(z + l)...(z + k-l). 

Replacing z by z + hq and applying Theorem [21 we find that the left hand side of J2J 
vanishes for $R(z) > no, hence by the identity theorem it vanishes in C slit along s rays 
from Ui to -co. We thus have 

d s 

= p (z) + (z- ...(*- u s ) a ° JT ((z ~u l + 1)~Y* 

k=l i=l 

= :;Po(z) + (z- Ul ) Ql G(z). 

If we assume w.l.o.g. that the Ui are in descending order w.r.t. the size of their real 
parts, then G{z) is analytic at z — u\. Since po(z) is not identically zero, this implies 
that (z — ui) ai is meromorphic at z = Ui, hence a± is an integer, a contradiction. The 
converse implication is trivial. □ 
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Another sequence that can be treated in a similar fashion is n an . So far, its non- 
holonomicity was only known for rational a 

Proposition 9. For every nonzero complex a, the sequence (n an ) n >x is not holo- 
nomic. 

Proof. Let F(z) be as usual. It is easily verified that (z + l) -Q ( z+1 ) F(z + 1) satisfies the 
growth conditions of Theorem|2| Appealing to condition (Al) completes the proof. □ 

Now let us consider rational functions in n and e™ with real coefficients. By Kho- 
vanskh's theorem and some commutative algebra, we can precisely characterize the 
holonomic sequences of this kind, even if they contain several exponentials. 

Theorem 10. Let atx, ■ ■ ■ ,a m be positive real numbers such that logai, . . . , loga m 
are linearly independent over Q. If f £ M(x, af , . . . , a^J is such that the sequence 
(/(n)) n >o is holonomic, then f(x) S M(x)[af, a± x , . . . , a^, <x^\. 

Proof. We write y = (y\ , . . . , y m ) for a vector of indeterminates and take some element 
g(x,y) — A(x,y) / B(x,y) of K(x,y) with the property that there is no such that 
(g(n, a", . . . , a^)) n > rao is a holonomic sequence. W.l.o.g. the polynomials A and B 
are coprime, and B depends on y m . By the closure of holonomic sequences under 
multiplication, we may also assume that B is irreducible. Moreover, B ^ R[a;]UR[a:]y m , 
or else there is nothing to show. Since this implies that Ay m and B are coprime, too, 
there are polynomials L(x,y) and M(x, y) such that 

L(x, y)A(x, y)y m + M(x, y)B{x, y) = R(x, y 1} . . . , y m -x). (4) 

Here, R is the resultant of the two polynomials Ay rn and B, viewed as polynomials in 
the variable y m . Dividing by B, we may thus assume that the numerator of g does 
not depend on y m . Pulling out the coefficient of the highest power of y m , we may at 
last assume that our g is of the form 

R(x,yi,...,y m -i) 

h -\ \-h D -lVm +Vm 

where D > 0, hj £ K(x, yi, . . . , y m _i), ho ^ 0, and the denominator is irreducible. 
Suppose now that g(n, a", . . . , a™ ) satisfies for large n. By Corollary (in fact, 
its obvious extension to [no, oo[), it follows that 

d d 

= ^p k {x)g{x + k, < +fc , . . . , a" +fe ) = ^p k {x)g k {x, of, ... , a x m ) 

for large real x. Here, we have defined 

9k(x, y) := g(x + k, a^yi, a k m y m ). 

The multiplicative independence of the numbers ax,..., a m implies the algebraic in- 
dependence of the functions af , . . . , a^, hence 

d 

^2pk(x)g k (x,y) = 0. 

k=0 
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Since the numerators of the gk do not depend on y m , there must be k > such that the 
denominators of go and gk, as polynomials in R(.t, y±, . . . , y m -i)[2/m]> have a common 
root in the algebraic closure of M.(x, yi, . . . , y m —i)- (Essentially, we are once more 
appealing to condition (Al).) But as these denominators are irreducible, the only way 
that they can share a root is if they coincide, up to a multiplicative "constant" from 
R(a;, yi, . . . , y m -i)- The leading coefficients w.r.t. y rn of the denominators of go and 
gk are 1 and , respectively, hence identifying the coefficients of yields 

ho(x,yx, . . .,y m -i) = a m kD h (x,a'ly 1 , . . . , a£,_ x y m -i)< 

But this equation cannot hold, since comparing the coefficient of any monomial (with 
non-zero coefficient) gives rise to a multiplicative dependence relation between the 
on- □ 

Singer and van der Put |24| have shown that the reciprocal (l/a„) of a holonomic 
sequence (a n ) is not holonomic unless (a„) is a periodic interlacement of hypergeo- 
metric sequences. Thus, in the special case m = 1 we can close the above proof after 
showing that the numerator of g(x, j/i) is w.l.o.g. independent of af. 

To see that the linear independence assumption in Thcorcm llOl is necessary, suppose 
that 

ci logai H h c m loga m = 0, 

where c\, . . . , c m are integers, and the are w.l.o.g. ordered such that 

ci, • • ■ , c fc > and c k +i, ■ ■ ■ , c m < 0. 
Then the sequence g((n, a™, . . . , a")) n >i corresponding to 

1 



1 i Ci Ci- — Ck+l —C m 

i + y x • • ■ y fc * - y fc+ i • ■ ■ 

is the constant sequence (l)„>i. 

5 Miscellaneous Results: Derivatives, Algebraic 
Functions, and Interlacing 

In this section we collect several results, mostly consequences of Khovanskh's theorem. 

Take some functions f(x) and g(x) defined for x > such that (f(n) n >i and 
(g(n) n >i are holonomic sequences. Which operations preserve holonomicity? Clearly, 
the sequences arising from f(x)+g(x) and f(x)g(x) are holonomic; we have been using 
this fact all the time. Since we assume that f(x) is defined for all positive real numbers, 
there are some other natural operations that would not make sense for sequences (a n ) 
in general. Integration does not necessarily lead to a function defining a holonomic 
sequence, as shown by f[x) = 1/x. Concerning differentiation, we note the following 
result. 

Proposition 11. Let f{x) be an elementary junction that is a quotient of D- finite 
functions defined for x > 0. If the sequence (f(ri)) n >i is holonomic, then the sequence 
(f'( n ))n>i is holonomic, too. 
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Proof. Assume that the sequence a n = f(n) satisfies ||TJ. By Corollary this implies 

d-l 

f(x + d) = -paix)- 1 Y,Pk{x)f{x + k), x> 1. 

k=0 

From this relation and Leibniz's rule it follows easily that, for every £ > d, the function 
f(x + £) is in the vector space spanned by f(x), . . . , f(x + d—l) and f'(x), . . . , f'(x + 
d — 1) over Since this vector space has dimension at most 2d, there must be a 

linear relation among f'(x), . . . , f'(x + 2d). □ 

Another natural question, which goes slightly beyond closed-form sequences, is 
whether sequences that satisfy some algebraic equation can be holonomic. It is clear 
that periodic interlacements of rational functions are holonomic, but no other algebraic 
holonomic sequences are known. For instance, we know from Section [3] that the se- 
quences n a with a £ Q\Z are not-holonomic, a fact already established by Gerhold [Hj 
with a number-theoretic argument. 

Theorem 12. Let f(z) be an algebraic function that is analytic in a neighborhood 
o/[l,oo[. Then the sequence (f(n)) n >i is holonomic if and only if f{z) is a rational 
function. 

Proof. Let F(z) = F(z; f,po, ■ ■ ■ ,Pd) be as usual. Let V be the algebraic variety 
generated by {(z,F(z)) : z > 1}. By Bezout's theorem, the points (z,F(z)) stay 
in one irreducible component of V for large real z (otherwise two components would 
have infinitely many intersections). As we suppose that F{z) vanishes for z £ N, this 
component has infinitely many intersections with the real axis, hence it equals the 
real axis, again by Bezout's theorem. But this shows that F(z) vanishes identically. 
Now observe that f[z) has an analytic continuation to a slit plane with finitely many 
singularities and branch points. If all the singularities are poles, then f(z) is rational 
p. 218], and we are done. Otherwise, we choose a singularity zq, not a pole, such that 
zq + 1, zq + 2, . . . are either poles or not singularities of f(z). Then we can appeal to 
condition (Al), since F(z) is not meromorphic at z = zq. □ 

A weaker assumption would be that the points (n, a n ) satisfy some algebraic equa- 
tion. The result should then be that a holonomic sequence arises only for a periodic in- 
terlacement of rational functions. The problem with the above proof is that the points 
(n,a„) may lie in several components of the variety generated by {(n,a n ) : n S N}. 
We can show as above that each of these components gives rise to a rational function, 
but what is missing is an argument that these rational functions must be interlaced in 
a periodic way. The special case where the rational functions are constants is doable, 
though: 

Proposition 13. // a holonomic sequence has only finitely many distinct values, then 
it is eventually periodic. 

Proof. It is known ^3 Exercise VIII. 158] that this holds if the sequence satisfies a 
linear recurrence with constant coefficients. We show that the result for a sequence 
{a n )n>i satisfying a recurrence of the form Q follows from this special case. Let N 
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be the maximum of the degrees of the p^s, and denote the coefficient of ri* in Pk{n) 
by bkj- Then 

d N 

= ^2p k {n)a n+k = y^£ 3 (n)n J 

fc=0 j=0 

for n > 1, where 

d 
fe=0 

By assumption, the set of N + 1 dimensional vectors 

£ = {(L {n),...,L N (n)) : n > 1} 

is finite. Since a non-zero polynomial has only finitely many roots, each non-zero 
vector from C can occur only for finitely many n. In particular, Ln{tl) must vanish for 
large n. By definition of N, some bjjy is non-zero, which yields a non-trivial constant 
coefficient recurrence for (a n ). □ 

If we assume in addition that the values a n in the preceding proposition are integers, 
then it follows from the Polya-Carlson theorem 2 that (a n ) must be a linear recurrence 
sequence. This observation is due to R.P. Stanley [2*2~], 

Finally, we briefly address the question whether non-holonomic sequences can be 
put together to produce a holonomic sequence. 

Proposition 14. Suppose that the sequence a n is a (not necessarily periodic) inter- 
lacement of the sequences logn and ^fn. Then a n is non-holonomic. 

Proof. Suppose that a n satisfies a recurrence of order d. Let us color the indices n 
where log n occurs black and the others white. In this coloring of N, there must be some 
pattern of length d + 1 that occurs infinitely often. Hence one of the 2 d+1 functions 

has infinitely many zeros. By Corollary |21 or Corollary it vanishes identically but 
this can once again be refuted by condition (Al). Hence a n is not holonomic. □ 



6 The Values of the Riemann Zeta Function at Odd 
Integers 

In this section we consider the sequence (C( n ))n>2, where £(s) = 53/>i denotes 
the Riemann zeta function. None of our conditions for property (A) is applicable. In 
particular, the zeta function is meromorphic in C, violating condition (Al). We are 
thus prompted to look for other arguments. If (C(^))n>2 was holonomic, then so would 
be the subsequence (£(2n)) n >i, which has the well-known representation 

C(2n) = ( _ ir «^J^ 
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in terms of the Bernoulli numbers, whose exponential generating function 



(2n)\ e z -l 

n>0 v ' 

is the reciprocal of a holonomic function. Holonomic functions are not closed under di- 
vision, and in fact Harris and Sibuya [S] have shown that the reciprocal 1/ f(z) of a holo- 
nomic function f(z) is holonomic if and only if the logarithmic derivative f'(z)/f(z) 
is algebraic. From this we infer that the Bernoulli numbers, hence (£(2n)) n >!, and 
hence (C( n ))n>2, are not holonomic. 

This straightforward argument tells us nothing about the values of the zeta function 
at odd integers. These numbers have attracted the attention of number theorists for 
(at least) the last decades, the most well-known result being Apery's proof [231 
of the irrationality of C(3). We will now show that the sequence (C(2ra + l)) n >i is 
not holonomic. Considering the series for £(2n + 1) summand by summand, we will 
deduce a full rank set of linear relations for the pk(n). As an analogy with the theory of 
transcendental numbers, we mention that the proof is evocative of the well-known fact 
that a rapidly converging series can be very useful for establishing the transcendence 
of a number. 

Proposition 15. The sequence (C(2n + l))„>i is not holonomic. 
Proof. Suppose that there are polynomials Pk(n) such that 
d 

^p fc (n)C(2n + 2fc + l) = 0, n > 1. 

By the estimate 



e>L+i 

with arbitrary e > 0, this implies 



L 



$> fc (n) = - J>(n) £ r 2 "" 2 *- 1 = 0((2 - e)- 2 "), (5) 

fc=0 k=0 l>1 

hence the polynomial on the left hand side of (JSJ is zero. We iterate this argument as 
follows: From 

d d 

5>fc(r»)(l + 2- 2 "- 2fe - 1 ) = - 5> fc (n) J2r 2n - 2k ~ 1 = 0((3 - e)^ 2 ") 

fc=0 k=0 l>3 

we obtain 

E^ = 0, n>l, 

fc=0 

and, continuing inductively, we find that 

d 



2^-^ = °' i,n>l. 



k=0 
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Since the Vandermonde matrix (£ 2k )i<e<d+i,o<k<d is invertible, the polynomials 
Po(n), . . . ,Pd{n) vanish at all n > 1, hence they vanish identically. □ 

7 Comments and Remarks 

When trying to establish property (A), we have to show that a certain function F 
does not vanish identically. A possible limit to this endeavor is hinted at by the 
celebrated results of Richardson J7| and Caviness [3]. Recall that they have shown 
that it is an algorithmically undecidable problem to determine whether a given function 
from certain classes of elementary functions vanishes identically. On the one hand, 
their functions have the absolute value function as a building block, which we do 
not consider; on the other hand, they assume a closed form expression, while our F 
contains arbitrary polynomials. 

As for property (B), we remark that there is an extension of Carlson's classical 
result due to Malliavin and Rubel [HI HOI • 

Applying Khovanskii's theorem only for quotients of D-fmite functions severely 
hampers the power of the theorem. Maybe the order of zeros of an elementary func- 
tions is bounded, which would make this restriction superfluous. An observation that 
might prove useful is that all elementary functions are differentially algebraic, i.e., 
P(z, /,/',..., f^) — for a nonzero polynomial P (in k + 2 variables). Hopefully 
this representation can be used to get an upper bound on the order of a zero of /. 
However, examples such as (/') 2 +/•/" = show that for algebraic differential equa- 
tions the boundedness cannot be established by the same straightforward argument as 
for D-finite functions. Rubel ^1 Problem 20] Problem 29] asks which sequences 
(zx, #2 1 • • ■ ) (counting multiplicities) can serve as zero sets of entire differentially alge- 
braic functions and remarks that the case (1, 2, 3, ... ) is still undecided. 

Furthermore, we remark that, surprisingly, in the case of integer sequences our 
restriction to sequences defined via elementary functions is not that restrictive after 
all. Laczkovich and Ruzsa |12j prove that every integer sequence can be represented 
in the form (f(n)) n >i, where / is a so-called naive-elementary function. However, this 
representation relies on the infinitude of zeros of the sine function, conflicting with the 
requirements of Khovanskii's theorem. 

We conclude our final remarks with four open problems. 

(i) Besides those described in Theorem^] which other elementary (in Khovanskii's 
sense) functions define holonomic sequences? (Note that (sinan) n >o clearly is 
holonomic.) 

(ii) Is a sequence satisfying an algebraic equation holonomic if and only if it is a 
periodic interlacement of rational functions? 

(iii) Can our method be pushed forward to accomodate recurrences of more general 
shapes and multivariate sequences? 

(iv) Let f(x) be a non-rational elementary function such that (f(n))n>i is holonomic. 
If s(x) is some function such that (f(s(n))) n >i is also holonomic, does it follow 
that s(x) is linear? 
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